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Abstract 

We study the energy flow between a one dimensional oscillator and a chaotic 
system with two degrees of freedom in the weak coupling limit. The oscillator's 
observables are averaged over an initially micro canonical ensemble of trajectories 
of the chaotic system, which plays the role of an environment for the oscillator. 
We show numerically that the oscillator's average energy exhibits irreversible dy- 
namics and 'thermal' equilibrium at long times. We use linear response theory to 
describe the dynamics at short times and we derive a condition for the absorption 
or dissipation of energy by the oscillator from the chaotic system. The equilibrium 
properties at long times, including the average equilibrium energies and the energy 
distributions, are explained with the help of statistical arguments. We also check 
that the concept of temperature defined in terms of the 'volume entropy' agrees 
very well with these energy distributions. 
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1 Introduction 



Low dimensional chaotic systems can, under appropriate circumstances, play 
the role of thermodynamical heat baths [1,2,3,4,5,6,7]. If a slow system with 
few degrees of freedom is weakly coupled to a fast chaotic system, the slow 
system's average trajectory can dissipate energy into the chaotic one at short 
times. 

One of the initial motivations for the consideration of low dimensional chaotic 
dynamics as environments for macroscopic systems was the work by Brown, 
Ott and Grebogi [8] on the ergodic adiabatic invariant. For a chaotic Hamilto- 
nian system with a slowly varying parameter, the volume of the energy shell 
is the only adiabatic invariant. Brown et al showed that the first order correc- 
tion to this invariant has a diffusive temporal behavior. Besides, if such time 
varying parameter is thought of as a second system coupled to the chaotic 
one, then this diffusive correction of the adiabatic invariant would lead to a 
dissipative force on the slow system [1] . 

This problem was again reformulated by Berry and Robbins [2] in terms of 
a system of interest interacting with an environment, or 'thermal bath'. The 
average force acting on the system of interest due to the coupling with the 
chaotic system can be calculated in adiabatic approximation. The lowest order 
part of this force was shown to be the Born-Oppenheimer force and the next 
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order to be a geometric magnetism type of force plus a deterministic friction - a 
force proportional to the slow system's velocity. Friction is therefore generated 
in the context of small systems, without the need for the thermodynamical 
limit, and the chaotic nature of the motion becomes the essential ingredient. 
Indeed, typical correlation functions of chaotic systems decays exponentially 
with the time, as opposed to the quasi periodic behavior observed in regular 
systems. Jarzynski [4] showed later that coupling to a low-dimensional chaotic 
motion can also lead to a 'thermalization' of the system of interest, very much 
like the thermalization of a Brownian particle interacting with a large thermal 
bath. Finally, Carvalho and de Aguiar [3] estabhshed a connection between 
de formalism developed by Caldeira and Leggett [9] for describing dissipation 
via coupling with a thermal bath and via coupling with a chaotic system. 

In this paper we revisit this problem from the classical point of view and 
consider specific examples of two-dimensional chaotic systems coupled to a 
one-dimensional oscillator. Our main purpose here is to understand the en- 
ergy fiow between the system of interest, which we shall call 'the oscillator' 
and the chaotic system to which this oscillator is coupled. We study the energy 
flow at short times and the approach to equilibrium at long times. In order to 
consider the chaotic system as playing the role of an environment, we assume 
that the only information available from this system is its initial energy. For 
the oscillator this implies that only microcanonical averages of its observables 
(over the chaotic system variables) are accessible. Therefore, a typical numer- 
ical calculation corresponds to fix an initial condition for the oscillator and 
to evolve an ensemble of trajectories whose initial conditions for the chaotic 
variables are randomly selected at a fixed (chaotic system) energy shell. 

For short times our numerical simulations show that the average energy of 
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the oscillator may increase or decrease, absorbing energy from the chaotic 
system or dissipating energy into the chaotic system. The initial energies of 
both systems is what ultimately dictates which of the two possibilities actu- 
ally occurs. In particular, there exists initial values of these energies such that 
no exchange occurs on the average. We use Linear Response Theory to study 
the energy flow in the short time limit. We show that the average motion of 
the oscillator follows a Langevin type of equation with frequency-dependent 
friction and a quadratic correction to the oscillator potential, similar to the 
Born-Oppenheimer force that appears in the adiabatic theory. We also derive 
a simple condition for dissipation or absorption of energy by the oscillator 
involving the ratio of the initial energies of the systems. This theoretical pre- 
diction agrees well with our numerical calculations for short times, but it is 
not accurate to predict the long time behavior. 

Our simulations show that, at long times, the average energy of both the os- 
cillator and the chaotic system tend to an equilibrium. The value of these 
equilibrium energies depend once again on the initial conditions. The connec- 
tion between asymptotic thermalization and initial conditions is well known 
for a Brownian particle. In that case, the increase or decrease of the average 
energy of the particle depends on its initial energy Eq and on the tempera- 
ture T of the thermal bath. The particle absorbs energy from the reservoir 
if £"0 < ksT and loses energy into it if £^0 > ksT, thermalizing always at 
ksT. Here we have a similar situation, with the increasing or decreasing of 
the average energy of the oscillator depending only on its initial energy and 
on the initial energy of the chaotic system. However, contrary to the case of 
the Brownian particle, the condition for equilibrium at long times is generally 
different from the condition of no energy exchange at short times. Despite 



4 



the theoretical work of Jarzynski [4] on the long term thermalization of these 
systems, it is still not clear how the asymptotic states depend on the initial 
conditions of both sub-systems. The energy distribution of the sub-systems 
at equilibrium is also an important open issue. These distributions are not 
likely to be of the Boltzmann type, since the systems are small, and they may 
depend not only on the initial conditions but also on the density of states of 
the systems involved [10,11,12]. In this paper we shall derive the long time 
equilibrium conditions and energy distributions explicitly for two model sys- 
tems. Finally, we use the definition of temperature proposed in [13] for small 
systems and check that it agrees completely with the statistical calculation in 
terms of the density of states and energy distributions. 

We emphasize that our approach uses the Linear Response Theory, and no 
explicit assumptions on adiabatic properties of the oscillator is required. De- 
spite the difficulties involved in the calculation of the response function for 
microcanonical ensembles, the formulation of this problem in terms of Lin- 
ear Response Theory is of great interest for the study of its quantum analog 
[7]. We recall that the usual formulation of quantum dissipation [9] involves 
response functions and that the adiabatic approximation leads to frustrating 
results in the quantum formulation [2] . The present work has some similarities 
with that of ref.[3], which also considered the model of an oscillator coupled 
to a chaotic bath to study dissipation at short times. Here we study both the 
short and long time limits, showing that the coupling with the chaotic system 
may indeed lead to an equilibration very much like that caused by the coupling 
with a large thermal reservoir. 

The outhne of the paper is as follows: in section II we describe our model 
systems and compute the correlation functions that are relevant for the calcu- 
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lation of the response functions. In Sec. Ill we review the classical theory of 
linear response and in Sec. IV we apply the theory to our model. We calculate 
the response functions explicitly in the microcanonical ensemble. In Sec. V we 
compute the average energy and average equation of motion of the oscillator. 
We also derive a condition for the initial increase or decrease of the average 
energy of the oscillator for short times, and compare it with numerical results. 
In Sec. VI we calculate the long time equilibrium values of the average energies 
using statistical arguments and compare them with numerical values. Finally 
in section VII we summarize our conclusions. 

2 The Model and First Numerical Results 

Our model Hamiltonian consists of a one-dimensional harmonic oscillator 
Ho{z) coupled to a chaotic system with two degrees of freedom Hc{x, y): 





where 




(2) 



and 



Vi{x, z) — ^xz. 



(3) 



For the chaotic Hamiltonian we consider two systems: 



Hc{x, y) 



(4) 
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and 
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(5) 



The Hamiltonian (4) is known as the Nelson system (NS) [14] and we shall refer 
to (5) as the Quartic system (QS) [20]. The NS exhibits soft chaos and is fairly 
regular for £'c^ 0.05, strongly chaotic for £'c^ 0.3 and mixed for intermediate 
values of the energy. The QS is invariant under a scaling transformation of the 
coordinates and the time, which implies that the dynamics on all energy shells 
are equivalent to each other. We shall explore this property later. The QS is 
integrable for a — 1.0, strongly chaotic for a ^0.1 and mixed for intermediate 
values of the parameter. 

We want to investigate the situation in which the chaotic system plays the 
role of an external environment for the oscillator. Therefore, we assume that 
detailed information about the chaotic system is not available. If the environ- 
ment were modelled by a heat bath, the only macroscopic relevant information 
would be its temperature. In the present case we assume that the only informa- 
tion available is the initial energy of the chaotic system. For the oscillator this 
implies that only microcanonical averages of its observables (over the chaotic 
system variables) are accessible. 

In order to characterize the two chaotic systems Eqs.(4) and (5), we show in 
Fig. 1 the NS correlation functions {x{0)x{t))e and {px{0)x{t))e, for Ec — 0.38 
and, in Fig. 2, the QS correlation functions for a — 0.1 and Ec = 5.0. These 
functions play important roles in the linear response theory to be developed 
later. The brackets (.)e stand for the average on the microcanonical ensemble 
of initial conditions. The correlation functions are obtained numerically with 
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a fixed time step symplectic integration algorithm [17] applied to the isolated 
chaotic system. In our calculations we used time steps of the order of 10~^ for 
figures 1, 2, 5, 6 and 7 (short times) and 10~^ for figures 3 and 4. The number 
of initial conditions in each case is indicated on the figure captions. 

These correlation functions can be well fitted by the expressions 

{x{0)x{t))e = cre~"* COS a;t, 

{Px{0)x{t))e = //e"^* sin Ot, (6) 

(see Figs.l and 2), which have the proper parity of the corresponding correla- 
tion functions. For the NS the fitting furnishes the decay rates a — 0.0418 and 
P = 0.0456, the amplitudes a = 1.865 and /j, — 0.409 and the frequencies of 
oscillation u = 0.1963 and = 0.2043 with ^ 10"^. For the QS we obtain 
a = 0.207, f3 = 0.208 and A = 0.206; a = 2.268, n = 3.67 and k = 4.10; and 
u = 1.1027, Q = 1.1481 and i/ = 1.189 with ~ 10"^- This type of correla- 
tion function is typical of chaotic systems in general [2,16] and also occurs in 
non-isolated systems subjected to random noise, such as a Brownian particle. 
The source of the correlation loss is, however, different in each case. Notice 
that the numerically calculated exponents a, (5 and A are all very similar and 
so are the frequencies a;, Q and v. 

When the coupling between the chaotic system and the oscillator is turned on, 
the overall conserved energy flows from one system to the other. The oscilla- 
tor's energy, in particular, fluctuates as a function of the time for each speciflc 

trajectory. The oscillator's average energy is calculated by taking an enscm- 
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ble of initial conditions uniformly distributed over the chaotic energy surface 
Ec- For the oscillator we fix only one initial condition, which we choose to be 
z{0) = and Pz{0) = 'j2mEjQ). The microcanonical ensemble of initial con- 
ditions is propagated and, at each instant, is calculated for each trajectory 
and its average value is computed. Figs. 3 and 4 show the oscillator's average 
energy {Eo{t)) as a function of the time for different values of Eo{0). We ob- 
serve that the average effect of the interaction with the chaotic system leads 
to a 'thermalization' of {Eo{t)). We also see that Eo{0) plays an important 
role on the long term behavior of {Eo{t)). As we shall see later using linear 
response theory, the ratio Eo{Q) / Ec{0) defines the initial rate of energy varia- 
tion in time. Zooming the short time behavior of the curves in Figs. 3 and 4 
(see Figs. 5 and 6) shows fast oscillations of {Eo{t)). As we shall see, part of 
these oscillations are due to a change in the effective potential acting of the 
oscillator, similar to the Born-Oppenheimer force of the adiabatic theory (see 
eq.(48) in section 5). Notice that the oscillator's equihbrium energy is different 
for each situation, which clearly distinguishes the small chaotic environment 
from an infinite thermal bath. 

As a last remark we note that, since the environment has few degrees of 
freedom, the oscillator's motion for a single realization (single initial condition) 
exhibits large fluctuations with respect to the average. 

3 Linear response theory 

The calculation of averages such as {Eo{t)) involve the calculation of the mi- 
crocanonical distribution function p{q,p;t) whose initial value is p{q,p;0) — 
6{H,{q,p) - Ee(0))/S(E,(0)), with S(E,(0)) = Jdqdp6{H,{q,p) - EM)- ^ 



9 



the chaotic system were isolated, p would be an invariant distribution and 
p{q,p; t) = p{q,p; 0). The coupling, however, causes the value of Hc{q{t),p{t)) 
to fluctuate in time, distorting the energy surface He = Ec{0) [15]. Linear 
response theory provides a way to calculate the first order corrections to this 
distribution in the limit of weak coupling [18,19]. 

Consider a Hamiltonian H{q,p) perturbed by a term of the form 

H, = -Aiq,p)X{t), (7) 

where A is an arbitrary function of the coordinates and momenta and X{t) 
is a function of the time. A generic distribution function p{q,p) follows the 
Liouville equation 

dp{t) 



dt 



i[L + L,{tMt), (8) 



where the Liouville operators are given by 



dHdp _ mdfA _ 
dq dp dp dq I 



(9) 



where the summation is taken over the whole set of canonical variables. 
The differential equation (8) has an integral representation of the form 



p(t) = e^(*-*°)^p(io) + J dse^(*-^)^iLi(s)p(s). (10) 

to 

Even if p(to) is an equilibrium distribution of H (canonical or microcanonical, 
for example), it will not remain invariant if the perturbation Hi is present. 
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However, if the perturbation is small, we may expand Eq.(lO) to first order in 
iLi to obtain 



p{q,p;t) = pe(g,p) + / dse^(*-^)^{i/i(s),pe(g,p)}, 



(11) 



where pe denotes the initial equilibrium distribution and p{q, p; t) is the time 
dependent non-equilibrium distribution up to first order in the perturbation 
Hi{t). Within this approximation, the average value of a general function 
B{q,p) can be computed as 

t 

{B)(t) = I dqdppe(q,p)B(q,p) + J ds(l>BA(t - s)X(s), (12) 

to 

where the response function (t>BA{t) is given by [18] 



with q{t) — e'*^g and p{t) — e**^p. Finally, expanding the Poisson brackets 

and integrating by parts we obtain 



The application of this theory to our model is as follows: we consider Vi 
as a perturbation on ifc, the coordinate z va. Vi being the analog of X{t) 
in (7). The initial microcanonical distribution for the chaotic system is then 
propagated according to Eq.(ll). Averages of observables involving the chaotic 
system variables can also be calculated with the help of Eq.(12). In addition 
to the consideration that the oscillator acts as a perturbation to the chaotic 



4>BA{t) = j dqdpB{q,p)e'''^{-A{q,p),p,{q,p)} 
= J dqdp{peiq,p),A{q,p)}B{q{t),p{t)), 



(13) 



<f>BA{t)^ I dqdppe{q,p){A{q,p),B{q{t),p{t))} 
^{{A{q,p),B{q{t),p{t))})e. 



(14) 
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system, we also consider the response of the chaotic system as a perturbation 
to the oscillator. This allows for the perturbative calculation of the oscillator's 
average trajectory and energy as well. 



4 Application to the model 



According to Eqs.(l)-(3), the equation of motion for the oscillator is given by 

z{t) + u;lz{t) = -—x{t) . (15) 
m 



This equation can be solved explicitly in the form 

t 

z{t) = Zdit) - — / dsT{t - s)x{s) (16) 



where Zd is the decoupled solution and T{t — s) = sin [a;o(t — s)]/<^o- 
Differentiating (16) with respect to t and multiplying by m we obtain 

t 

Pz{t) = Pzait) - 7 / dsxit - s)x{s) (17) 



with x(t — s) = cos [u!Q{t — s)]. Taking z{0) = the decoupled solutions Zd{t) 
and Pz^it) become 

Zd(t) = ^ sin (uot) p,,(t) = p,(0) cos (uot) (18) 
mojQ 

with ^2(0) — Y^2m£'o(0). The oscillator's average energy {Eo{t)) can be com- 
puted as 

^EM)-^^^^^^ (19) 

12 



with 



{Pl(t)) = pl(t) - 27p,,(i) I dsx(t - s){x(s)) 



+ 



t t 

,,2 



rjdsj duxit - s)x{t - u){x{s)x{u)) (20) 





and 



(z'it)) = zlit) - ^z,{t) j dsr{t - s){xis)) + 

TJX J 


2 * * 

^ I ds I duVit - s)V{t - u)(xis)x{u)). (21) 







Notice that the quantity {Eo{t)) is very different from Eo{{z), ipz)), which is 
the energy of the average trajectory, considered in [3]. 

Linear response theory now furnishes 

t 

{x(t)) = {x(t))e - 7/ ds(l)^x{t - s)z{s) (22) 



where (f)xx{t) — {{x{0), x{t)})e is the response function given by Eq.(14) with 
A = B = X. Since Hc{—x) = Hc{x) for both chaotic systems given by Eqs.(4) 
and (5), {x{t))e — 0. Substituting (22) into (20) and (21) and considering 
terms up to order 7^ we obtain 

= pU^) + 27' lo dsx(t - s) du(l)xx{s - u)zd{u) 

(23) 

+7^ Jo Jo dux(t - s)x(t - u){x(s)x{u))e, 
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{z\t)) = zl{t) + 1^ f^dsT{t - s) S^du^Us - u)za{u) 



(24) 



+^ Jo* ds /o* duT{t - s)T{t - u) {x{s)x{u)),. 



Eqs. (23) and (24) represent the average behavior of the harmonic oscillator 
when interacting with a generic x-symmetric system through the perturbation 
Vi. They show that all the influence of the interacting system is contained in 
the functions {x{Q)x{t))e and (pxxit)- If the interacting system is integrable, 
these functions exhibit quasi-periodic oscillations. For chaotic systems, on the 
other hand, they typically decay exponentially with time, leading to qualita- 
tively distinct average results. 

The response function is given by 4>xx{t) — ({^(0), a;(t)})e where {.} is the 
Poisson bracket with respect to the initial conditions [18] . Since Px{0), Py{0), 
x{0) and y{0) are independent variables we integrate by parts and obtain 




(25) 



where dV = dx{0)dy{0)dpx{0)dpy{0), 



Pe = 5(//e-^c(0))/S(Ee(0)) 



(26) 



is the microcanonical distribution and 




(27) 



is a normalization factor. We obtain 




E(E,(0)) 



1 



jdV 



OH, 



■Px{0)x{t). 



(28) 
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This integral can be performed explicitly by changing to variables across the 
energy shell, i^c, and along the energy shell, ^, 6 and ip: 



(29) 



where = {C,,9,ip), J = J{Hc,^,0,(p) is the Jacobian of the transformation 
and the derivatives are computed at Ec{0). Similar results for response func- 
tions calculated on the microcanonical ensemble can found in [21]. For the NS, 
the explicit transformation is given by 



/ , \l/2 2 

^"(otsj V -f^c COS ^ COS ^ y = ^ + V -f^c sin ^ COS </? 

(30) 

Px — \/2Hc cos ^ sin 9 py — ^/2H^ sin ^ sin (f 



with J = HcfN{Q)- The response function simplifies to: 

^"W-^^ + (al;W0M.))) . (31) 



For He close to Ec{0) we assume that the amplitude of the correlation function 
{px{0)x{t))e varies linearly with He and we make the approximation 

We obtain 

'^"^^^^ = e|o) ^^^^ 
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which is similar to the canonical response function [18]. 



For QS, the explicit transformation is given by 



X 



2Hc 


cos 6 


cos 29 


. Vl+a 


2Hc 


cos 6 


cos 20 





+ 



sin 



n-a 



(34) 



Px = y/2Hc cos if cos ^ 



Py = \/2HcSm(p cos^ 



with < ^ < ^0 < 7''/4 and tan^o = y (1 ~ + «)• The jacobian is given 
J = Hl^"^ fqiQ) and the response function (29) simplifies to: 



(35) 



The last term of (35) can be calculated exphcitly because the QS has an energy 
scaling dynamics [20] given by: 



Px(t)=(|^) pL(t'), ^W=(§) ^'(0, t=(§) (36) 



EL 



EL 



where p^, x' and t' represent the motion with the energy E'^ and Px^ x and t 
the motion with energy Ef.- 

Using (36) in (35), we obtain 



dE, 



3/4 



^E'J\0)H'J' 



P'AO)x'{t') + 



EM J 



3/4 



p'M 



dx'{t' 

"oh:. 



(37) 
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where the last term is given by 



Ec{0) J "^^^ ' dH^ \EM ) dt' dH^ 



Thus, 



(38) 



^ {pMm) + (39) 



and, from (35) and (39), 

The simphfication we obtained, due to the specific form of the Jacobian, is 
vahd for a large class of Hamiltonians, including the billiards. 



5 Short time dynamics 

In order to derive some explicit formulas, we consider the expressions (6) for 
the correlation functions derived in section 2 by fitting the numerical results. 
For the NS we assume fl lu and /3 a; a. Using Eqs. (19), (23) and (24) we 
obtain the following result: 

{E,{t)) = EM + ^{B + At + fit) + g{t)), (41) 
m 

where is a constant, f{t) is an oscillatory function and g{t) is proportional 
to e~"*. The important result is the coefficient A 



■[(a;o-a;)2 + a2][(a;o + a;)2 + a2]- 
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For fixed oscillator frequency cuq and a given chaotic energy shell Ec{0) (and, 
consequently, for given a, fi, uj and a), the ratio Eo{0) / Ec{0) is the responsible 
for the average increase or decrease of {Eo{t)). The short time equilibrium in 
the energy flow is given by the condition A — 0, or 



Ec{0) AlMjJ 



uj'^+uj' + a') . (43) 



We now turn to the equation of motion of z{t) under the average effect of the 
interaction with the chaotic system. Prom Eqs.(15) and (22) 



„2 * 



{z{t))+^lm) = ^-{m = - f d«0..(i - s){z(s)). (44) 



Integrating by parts yields 
7'F(0) 



m) + Ho - ^— ^ {m + ^-UsF{t- s){z{s)) + ^-z{Q)F{t) 440) 
\ m m J m 



where 



~ " £;^(0)(a2 + a;2) {^ cos [uj{t - s)] + asm [uj{t - s)]}. (46) 



Eq.(45) shows that the interaction produces a harmonic correction to the 
original potential, a frequency-dependent dissipative term and an external 
force proportional to ^(0). The choice z{0) = simplifies (45) and turns it 
into an average Langevin equation. 

We note that the correction to the harmonic potential is repulsive and may 
qualitatively change the average motion of the oscillator if 7^ > rmjQ/F{0). 
The critical perturbation 7^ — mujQ/F{Qi) is related to the stability of the 
equilibrium points of the full three-dimensional system, which becomes unsta- 
ble for large 7, establishing a direction of unbounded motion. For the NS, for 
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example, the origin (x, y, z) — (0, 0, 0) becomes unstable at 7^ = Q.lm/ujQ (see 
Eq.(4)). The function F(0), on the other hand, can be computed numerically 
from Eq.(46) and results in approximately 10. This value can also be estimated 
by noticing that F{t) has dimensions of {x{Q)x{t)) / Ec, which, at i = 0, yields 
2{x{0Y) / Ec = 1/0.1. Therefore one must choose 7 < 7c in order to consider 
the coupling as a perturbation to the original system, avoiding the introduc- 
tion of these instabilities. We remark that similar corrections are observed in 
the theory of Caldeira and Leggett [9] and in the adiabatic calculations of 
Wilkinson-Berry- Robins [2]. 

Fig. 5 shows a comparison between the numerically calculated 'bare' oscillator 
energy {Eo{t)), where 



^'^'^ ^^^^ 



the 're-normahzed' oscillator energy {Eor{t)), where 



ii,!=|£!_l!,(oy, (48) 



and the expression (41) without the oscillating term f{t). We have chosen 7 
and m so that 6i;q— 7^F(0)/m > 0. We also chose loq so that g{t) decreases very 
fast, i.e., e-"/'^o a; lO""^. In this case only the linear and the oscillating terms 
in Eq.(41) are important. We have subtracted the oscillating part of Eq.(41) 
in figure 5 to highlight the linear increase or decrease in the average energy. In 
the time scale of Fig. 5, which corresponds to several periods of the decoupled 
oscillator, the linear behavior describes very well the numerical results. Fig. 
5(b) shows the equilibrium situation according to Eq.(43). 
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For the QS the correlation functions can also be simplified with the assumption 
that P ^ a ^ X (we shall keep for now u; ^ Vl ^ u). We get from (19), (23) 
and (24) 

{E,{t)) = EM + ^{B' + A't + fit) + g'{t)), (49) 
m 

where B' is a constant, f'{t) is an oscillatory function and g'{t) is a sum of 
terms proportionals to e~"* and te~"*. The coefficient A' is given by 



^ — O U. . , .\2 , „,21 U. . , , .\2 , „,21 ' y^^) 

where 



, _ [(^0 - + o^'] [(^0 + + o^'] .... 



and 



^ , // [{cuq - + g^] [{up + O)^ + g^] 

^ = ' + 5^1^ ^ 

[(cug + ft:^)^ + 2;/^(cuo + /t)(cJo - k) - 3u^] 

[{uJo - vy + K^f [(cuo + i^)' + ^^^f 



(52) 



We note that if we set k = (which amounts to cancel the last term of (39)) 
and with a; = Q, (50) becomes very similar to (42). For k 7^ 0, can be 
positive or negative and if < 0, {EJ^)) cannot decrease with time. This 
possibility never happens for the NS. The equilibrium in the energy fiow is 
given by the condition: 

Eo% ^ 2Aa (a;2+a;2 + Q,2) 

£;e(o) 5/3Q 77 ■ ^ ' 

The equation (45) remains valid for QS, but the function F{t) is given by: 
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~4Ec(0)(A;2 + i/2) 



{kcos[u{t — s)] — i'sm[{u{t — s)]} 



^ {{u^ - k^) cos[v{t - s)] + 2vk sin[z/(i - s)]} 



+ 



4£;e(0)(A;2 + i/2)2 



4£;e(0)(Q2 + A;2) 



{f] cos[i/(t - s)] + A; sin[i/(t - s)]} . (54) 



Figure 6 shows the comparison between the numerically calculated {Eoit)) and 
{Eor{t)) and the expression (49) without f'{t). The parameters were chosen 
in such a way that e""/*^" ~ 10^^ and e^°'^^"/uJo ^ 10^'^, so that the linear 
and oscillating terms are the most important. Again we have subtracted f'{t) 
to highlight the linear behavior and we have chosen 7 and m so that cu^ — 
-f^F{Q)/m > 0. 

The agreement between numerical and analytical results for the equilibrium 
condition is only reasonable for the QS. The amplitude of the oscillations in 
{Eo{t)) are bigger than in the NS case because of the different energy scale 
and coupling strength {Ec{0) = 0.38, 7 = 0.006 for the NS and Ec{0) = 5.0, 
7 = 0.01 for the QS). In the approximation a — l3 — X and cu ^ D, ^ u, 
the equilibrium condition A' = gives Eo{0)/Ec{0) = 0.45, but if we use 
Lo — O, u, it changes to Eo{0)/Ec{0) — 0.55. This shows that small changes 
in the frequencies of the correlation functions change significantly the value of 
Eo{0) / Ec{0) for the equilibrium condition. Fig. 7 displays {Eo{t)) for others 
values of Eo{0)/Ec{0). It shows that the slope of the average energy as a 
function of the ratio Eo{0) / Ec{0) has a very shallow minimum, which justifies 
its poor determination with linear perturbation theory. Fig. 7 suggests that 
the equihbrium condition may actually be located between Eo{0) / Ec{0) — 2.5 
and Eo{0)/Ec{0) = 3.0. 
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6 Long time behavior and 'thermal' equiUbrium 

For long times the oscillator and the chaotic system reach an equilibrium, in 
the sense that their average energies tend to constant values. This is shown 
in Figs. 3 and 4 for the NS and the QS respectively. Notice that the energy 
ratio at the asymptotic equilibrium, at long times, is not necessary related to 
the equilibrium condition at short times, since non-linear effects are certainly 
important in the later. In order to characterize the equilibrium it is important 
to understand the energy distribution within each sub-system. It is particularly 
interesting to check if the oscillator follows Boltzmann's distribution. 

To obtain the distributions numerically, we construct histograms in which the 
values of and Ec are extracted from each trajectory of the ensemble for 
a fixed (long) time. The energy axis is divided into bins and the number of 
trajectories of the ensemble for which the oscillator's energy fall into each 
bin is counted. The same process is performed with respect to the chaotic 
system's energy. Figure 8 shows the energy distributions thus obtained for the 
oscillator and the chaotic systems, NS and QS. The vertical axis shows the 
counts (occurrences) in percents of the total number of trajectories. In both 
cases it is clear that the oscillator does not follow the Boltzmann exponential 
law. This is actually not surprising, since the chaotic system is small (with 
only two degrees of freedom) and thus its energy is comparable to that of 
the oscillator. The Boltzmann distribution comes out naturally only when the 
system of interest is in contact with a reservoir of much larger energy (spread 
among its many degrees of freedom) [23], a condition that is not fulfilled here. 
In this context, it is also not clear whether a temperature can be defined in 
the present situation. We shall return to this question in a moment. 
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The energy distribution of the subsystems can be totally understood in terms 
of their density of states. To see this, we assume that two hypothesis that 
are usual in the statistical physics of large systems [23] are also valid for our 
small systems. The first is to consider that, at equilibrium, all states of the 
full system, Eq.(l), are equally probable, i.e., the hypothesis of equal proba- 
bilities a priori. The second assumption concerns weak interactions between 
the subsystems. The number of states dN{E) of the full system, for which the 
oscillator has energy between E and E + dE, can be written in terms of the 
density of states 



where H[x,y,z) is given by Eq.(l) and dV = dxdpxdydpydzdpz- Neglecting 
the interaction potential Vj in Eq.(l), we can write dN{E) — no{E)nc{ET — 
E) dE, where rio and Uc are the oscillator and chaotic system densities of 
states, respectively, and Et is the total energy. 

Based on these two assumptions, we can calculate the probability that the 
oscillator has energy between E and E + dE. The probability density Po{E) is, 
because of the first assumption, proportional to the number of states dN{E): 




(55) 



Po{E)dE (X dN{E) 



no{E)nc{ET - E)dE. 



(56) 



Likewise, the probability density of the chaotic system is 



Pc{E)dE (X nc{E)no{ET - E)dE. 



(57) 



We find that no{E) is a constant and that nc{E) oc E for NS and nc{E) oc E'^^'^ 
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for QS. Thus we obtain 



Po{E) oc {Et - E) 



< 



for the NS 



(58) 



Pc{E) oc E 



and 



Po{E) oc {Et - Efl^ 



< 



for the QS 



(59) 



oc E^^ 



Of course these expressions are meaningful only for < E' < Et- 

The full lines in Fig. 8 show a linear fitting for the case of the NS and a square 
root fitting for the QS. The fittings agree very well with the numerical results. 
The sudden decrease of the numerical distributions of the chaotic systems for 
high energies is due to the constraint that Et is fixed. 

Finally, with Eqs.(58) and (59), we can also calculate the oscillator's average 
energy in the equilibrium and compare the results with the values in figures 3 
and 4. We have 



is Z = El/2 for the NS and Z = ^E^J^ /?, for the QS. We obtain E„ = Et/?> 
for NS and Eo = 2Et/5 for QS. From the probability densities Pc for N.S. and 




(60) 







where the normalization constant 




(61) 
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Q.S., we obtain Ec — 2Et/S and Ec — 3Et/5, respectively, for N.S. and Q.S.. 
For the parameters of Fig. 3a, ^^(0) = 0.38 and Eo{0)/Ec{0) = 1.0, we find 
Eo = 0.253 and for Fig. 3b, where ^^(0) = 0.38 and Eo{0)/Ec{0) = 0.1, we get 
Eg — 0.139. For Figs. 4a and 4b we obtain Eg — 12 and Eg — 2.4 respectively. 
In all cases the numerical values of the oscillator equilibrium energies is very 
close to the statistical prediction. 

We now return to the question about the temperature. We will show that not 
only it is possible to define a temperature for our subsystems but also that 
the equilibrium predicted by equating these temperatures results in the same 
partition of average energies predicted by the statistical analysis above. We 
first consider the usual definition of temperature given by 

]__dS_ 
f~dE 

where 

S = ke In n{E), (63) 

is the entropy and n{E) is the density of states given by (55). The thermal 
equilibrium between the oscillator and the chaotic system implies dSo/dEg — 
dSc/dEc- However, because no{E) does not depend on E, this gives = T^. = 
oo, and the equilibrium condition becomes useless. 

Recent studies [22,13] have proposed modifications in the calculation of the 
entropy that, although irrelevant for large systems, make significant differences 
for small systems. Ref.[22] suggests dynamical corrections to the Boltzmann 
principle (Eq.(63)). On the other hand, Ref.[13] argues that the entropy in 
Eqs.(62) and (63,) should be replaced by 
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(62) 



S^^kBln<5{E), (64) 

where n{E) — It has been shown that Eqs.(63) and (64) lead to identical 

results in the thermodynamic limit, but not for small systems, where (64) is 
still able to describe well the results of numerical simulations. 

For the present model we have ^o{E) oc E, ^ns{E) oc E'^ and ^qs{E) oc E^/"^. 
In this framework, the equilibrium condition can be obtained equating the 
temperatures of the oscillator and the chaotic system, computed with the 
modified entropy Eq.(64). For the NS we find 

dEo dEc ^ Ec 2 ^ ' 

and for the QS 

d\n^o{Eo) _ d\n^Qs{E,) E^ _ 2 



dEo dEc E^ 3' 



(66) 



These results are in complete agreement with those obtained via the proba- 
bility densities. Considering that our theoretical equilibrium energies describe 
very well the numerical calculations and the agreement between these energies 
and the thermal equilibrium conditions, we can conclude that the temperature 
is indeed a good parameter for characterizing the equilibrium. 

As a last remark we note that the initial distribution function, which is mi- 
crocanonical only in the chaotic degrees of freedom, is not expected to evolve 
to a fully microcanonical distribution over the entire system. This is because 
the dynamics of the full system is probably mixed, not ergodic. However, the 
agreement of the energy distributions at long times with the above calcula- 
tion suggests that the dynamics is at least 'close' to ergodic, in the sense that 
typical trajectories explore a large fraction of the available energy shell. 
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7 Conclusions 



Our numerical results show that the coupling of an oscillator to a low dimen- 
sional chaotic system simulates very well some aspects of a Brownian particle 
in a harmonic potential in the presence of a thermal bath. In particular, the 
average energy of the harmonic oscillator shows irreversible behavior and tends 
to an equilibrium at long times. For short time scales the average dynamics 
of the oscillator can be studied by linear response theory, where irreversibil- 
ity is seen to result from the temporal decay of the chaotic correlations. The 
average motion of the oscillator follows a Langevin type of equation whose 
frequency-dependent dissipation depends only the amplitude, frequency and 
decay rate of the chaotic correlations. 

The initial energies of the sub-systems play an important role in the prob- 
lem. Ec{0) defines the dynamical regime of the chaotic system, and the ratio 
Eo{0)/Ec{0) defines if the oscillator will absorb or dissipate energy to the 
chaotic system. This is very similar to the usual thermalization of a Brownian 
particle, with -E'c(O) playing the role of the average energy ksT of the ther- 
mal bath. An important difference, however, is that the average energy of the 
chaotic system is also affected by the coupling, which is a direct consequence 
of its small number of degrees of freedom. 

For long times the average energy of both oscillator and chaotic system tend 
to equilibrate. The value of the average energy and the energy distribution 
at equilibrium can be calculated assuming equal probabilities of the available 
microscopic states of the full system and weak coupling, so that probabilities 
over the full system can be approximated by the product of the probabilities 
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over each subsystem. 

It is interesting to see that both the short time and the long time analysis allow 
the determination of initial energies £"0(0) and Ec{0) such that no exchange 
of energies occur in the average. For short times it is given by the condi- 
tion A — 0, Eq.(43) for the NS system. For the parameters of Fig.3 it gives 
Eo{0) / Ec{0) ~ 0.25. For long times, since Eo = Et/S, imposing Eo = -E'o(O) 
and Et = Eo{0) + Ec{0), we find E„{0)/E^{0) = 0.5. The two estimates clearly 
disagree, which means that for intermediate times non-linear corrections to the 
linear theory become important and change the short time tendency of the 
average energy. In other words, the short time dynamics is completely de- 
termined by the properties of the isolated chaotic system, whereas the long 
time behavior is dictated by the statistical properties of the full system. The 
same reasoning apphes to the QS. Fig.9 shows {Eo{t)) for -E'c(O) = 0.38 and 
£"0(0) = £'c(0)/2 = 0.19. This is the condition for no exchange of energies at 
long times, but corresponds to a situation where the oscillator should dissi- 
pate at short times. And that is exactly what happens: the long time behavior 
displayed in Fig. 9 (a) shows that the oscillator's average energy is indeed ap- 
proximately equal to its starting energy. However, for short times. Figs. 9(a) 
and (b), it clearly dissipates energy, re-absorbing it back latter on. The dis- 
agreement between short and long time equilibrium conditions may also be 
related to the low dimensionality of the chaotic system [21]. 

The low dimensionahty is also reflected in the fact that the oscillator's equi- 
librium energy is different for each initial set up, clearly distinguishing the 
small chaotic environment from an infinite thermal bath. Besides, the oscil- 
lator's motion for a single realization (single initial condition) exhibits large 
fluctuations with respect to the average. 
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Finally, the temperature defined from the volume entropy [13] describes 
very well the equilibrium condition of our model, since it agrees with the 
equilibrium conditions found by the probability densities and, consequently, 
with the numerical results. We believe this result contributes to the question 
of '. . .whether the volume or the area entropy is the correct starting point of 
thermostatistics for small systems' [13]. However, it is important to remember 
that, although it is possible to define such a temperature, there are always 
large fluctuations around these equilibrium values because of the low number 
of degrees of freedom. Moreover, it is not clear why the probability densities 
are given by the density of states and the temperature by its integral. 

We note that the validity of the two assumptions we used here to describe 
the long time equilibrium regime are the subject of a debate that touches on 
the foundations of statistical mechanics [24]. The assumption of equal prob- 
abilities a priori is related to ergodicity and mixing. Here we have applied 
this condition without much justification, especially because the global sys- 
tem Eq.(l) is probably not mixing. However, the results we obtained from this 
assumption agrees very well with the numerical calculations. Thus, we believe 
that this type of low-dimensional models can be good testing grounds to the 
study of these topics. 
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Fig. 1. Correlation functions for the NS for Ec = 0.38: (a) {p^{0)x{t))e; (b) 
{x{0)x{t))e. The full line shows the numerical results and the dashed line shows 
the fitting. The averages were computed using 35000 initial conditions. 
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Fig. 2. Correlation functions for the QS for Ec = 5.0 and a = 0.1: (a) {px{0)x{t))e; 
(b) {x{0)x{t))e', and (c) and (j>x{0)Px{t))e- The full line shows the numerical results 
and the dashed line shows the fitting. The averages were computed using 30000 
initial conditions. 
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Fig. 3. Numerical results for {Eo{t)) with Ec{0) = 0.38 for NS: {si)Eo{0)/Ec{0) = 1.0 
and {h)Eo{0)/Ec{0) = 0.1. The oscillator's parameters were set to m = 200.0, 
uq = 0.005 and the coupling constant to 7 = 0.006. The oscillator's period is 
To ^ 1260. We used 62480 initial conditions. 
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Fig. 4. Numerical results for {Eoit)) with Ec{0) = 5.0 and a = 0.1 for the QS: 
{8:)Eo{0)/Ec{0) = 5.0 and {h)Eo{0)/Ec{0) = 0.2. The oscillator's parameters were 
set to m = 10.0, ujo = 0.01 and the coupling constant to 7 = 0.01. The oscillator's 
period is Tq 628. We used 40000 initial conditions. 
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Fig. 5. Average oscillator energy at short times with the NS as chaotic system. 
The dashed line shows {Eo{t)) and the doted line shows {Eor{t)), both obtained 
numerically The full line shows Eq.(41) without f{t). (a) Eo{0)/Ec{0) = 1.0, (b) 
Eo{0)/Ec{0) = 0.25 and (c) Eo{0)/Ec{0) = 0.1. The oscillator's parameters, cou- 
pling constant and number of initial conditions are the same as in Fig.3. 
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Fig. 6. Average oscillator energy at short times with the QS as chaotic system. 
The dashed line shows {Eo{t)) and the doted line shows {Eoj.{t)), both obtained 
numerically. The full line shows Eq.(49) without f'{t). (a) Eo{0)/Ec{0) = 5.0 and 
(b) Eo{0) / Ec{0) = 0.2. The oscillator's parameters, coupling constant and number 
of initial conditions are the same as in Fig.4. 
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Fig. 7. Average oscillator energy at short times with the QS as chaotic sys- 
tem. The dashed line shows {Eo{t)) and the doted line shows {Eor{t)), both ob- 
tained numerically. The full line corresponds to Eo{0). (a) Eo{0) / Ec{0) = 2.0, (b) 
Eo{0)/Ec{0) = 2.5, {c)Eo{0)/Ec{0) = 3.0 and Eo{0)/Ec{0) = 3.5. The oscillator's 
parameters, coupling constant and number of initial conditions are the same as in 
Fig.4. 
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Fig. 8. (a)Distribution of energies of the ensemble of 30000 trajectories at the time 
t = 8 X 10^, for the parameters of Fig. 3a. The full line shows to linear fitting, the 
dashed line shows Eg and dotted line shows E^, (b) Distribution of energies at the 
time i = 16 X 10^, for the parameters of Fig. 4b and 30000 trajectories. The full line 
corresponds to square root fitting, the dashed line to Eg and the dotted line to Ec- 
The energy is shown in units of Et/SO. 
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Fig. 9. (a) {Eo{t)) for EdO) = 0.38 and £^0(0) = Ec{0)/2 = 0.19. (b) Magnification 
of the time interval from to 2 x 10^ showing the initial dissipation of energy. 
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